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Abstract 

Activity studies of solar-type stars, especially with reference to the status of our current Sun among 
them, have exposed the importance of (1) homogeneously selecting the sample stars and (2) reliably eval- 
" uating their activities down to a considerably low level. Motivated by these requirements, we conducted 

an extensive study on the activities of 118 solar-analog stars (of sufficiently similar properties to each 
other) by measuring the emission strength at the core of Ca II 3933.663 line (K line) on the high-dispersion 
spectrogram obtained by Subaru/HDS, where special attention was paid to correctly detecting the chromo- 



spheric emission by removing the wing-fitted photosphcric profile calculated from the classical solar model 
atmosphere. This enabled us to detect low-level activities down to logi?' ~ —5.4 (R 1 is the ratio of the 
chromospheric core emission flux to the total bolometric flux), by which we could detect subtle activity 
differences which were indiscernible in previous studies. Regarding the Sun, we found logi?Q = —5.33 near 
to the low end of the distribution, which means that it belongs to the distinctly low activity group among 
solar analogs. This excludes the once-suggested possibility for the high frequency of Maunder-minimum 
stars showing appreciably lower activities than the minimum-Sun. 
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. 1. Introduction However, this result could not be confirmed by a simi- 

lar analysis done by Hall and Lockwood (2004), who re- 
It is of great interest for solar as well as stellar astro- ported based on many repeated observations of Ca II H 
physicists to compare the activity of our Sun to those of and K lines for 57 Sun-like stars along with the Sun at 
a number of other similar solar-type stars, since it may Lowell Observatory that such a bimodal distribution of 
provide us with an opportunity to infer the trend of solar S'-indcx (as suggesting the existence of a considerable frac- 
activity on a very long astronomical time scale (i.e., inves- tion of appreciably lower activity stars than the current 
tigating the long-time behavior of a star may be replaced Sun) is not observed; actually, even the S-values of 10 
by studying many similar stars at a given time) . "flat-activity" stars turned out to be comparable with (or 

Baliunas and Jastrow (1990) argued based on the re- somewhat larger than) the typical solar-minimum value, 
suits of Mt. Wilson Observatory's HK survey project for Furthermore, Wright (2004) pointed out an important 
74 solar-type stars that the distribution of 5-index (nearly problem in Baliunas and Jastrow's (1990) sample selec- 
equivalcnt to cx J line F\d\/ J .F^dX; i.e., the ratio of in- tion. He concluded by examining the absolute magni- 
tegrated core-flux of Ca n HK lines to the continuum flux; tudes of their sample based on Hipparcos parallaxes that 
cf. Vaughan et al. 1978) is bimodal, with about 1/3 show- many of those "Maunder-minimum stars" with consider- 
ing appreciably smaller activities than the Sun, which they ably low S indices are old stars evolved-off the main se- 
interpreted as being in the "Maunder-minimum" state of quence, which suggests that their apparently low activity 
activity. If this is true, it may mean that the spotless is nothing but due to the aging effect without any rele- 
phase of considerably low-activity such as occurred in late vance to the cyclic or irregular change of activity in solar- 
17th century in our Sun (Eddy 1976) may not necessarily type dwarfs. Thus, fairly speaking, the original claim by 
be an unusual phenomenon in the long run. Baliunas and Jastrow (1990) appears to be rather pre- 

mature and difficult to be justified from the viewpoint of 
these recent studies. 

Yet, the issue of clarifying the status of solar activity 
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among Sun-like stars does not seem to have been fully set- 
tled and further investigations still remain to be done: 

— First, since understanding the activity of the Sun from 
a comprehensive perspective is in question, comparison 
samples should comprise stars as closer to the Sun as pos- 
sible. Admittedly, those authors surely paid attention to 
this point: Baliunas and Jastrow's (1990) HK project tar- 
gets were in the B — V range of 0.60-0.76, while Hall and 
Lockwood's (2004) Sun-like stars sample were chosen from 
stars of 0.58 < B — V < 0.72, both narrowly encompassing 
the solar (B - V) Q of 0.65 (Cox 2000). However, the ho- 
mogeneity of these samples are not yet satisfactory. Could 
they be made up of further more Sun-like stars or solar 
analogs? 

— Second, it appears that the precision of detecting low- 
level activity has been insufficient. Although Mt. Wilson 
S index reflects the core-emission strength (equivalent 
width) of Ca n HK lines, it would not be a sensitive 
activity indicator any more, when the emission becomes 
weak, as it stabilizes at a constant value determined by the 
photospheric absorption profile. While an indicator for 
the pure-emission strength, i?^ K (= -Rhk — -Rphot! where 
-Rhk = Fn K /F hol and i? pho t = F phot /F ho {), has also been 
introduced to rectify this shortcoming and widely used, 
the photospheric component i? p hot is in most cases only 
roughly evaluated as a simple function of B — V (e.g., 
Noyes et al. 1984) and thus its accuracy is rather ques- 
tionable. Wright (2004) also pointed out the importance 
of correctly subtracting this component, in view of its pos- 
sible dependence on other stellar parameters (i.e., not only 
on B — V or T e g, but also on log g and [Fe/H]). 

These requirements motivated us to conduct a new 
investigation on this subject, since we have been en- 
gaged these years with the extensive study of 118 Sun- 
like stars selected by the criteria of 0.62 < B — V J$ 0.67 
and 4.5 < My < 5.1 (M ViQ = 4.82; Cox 2000), a quantita- 
tively as well as qualitatively ideal sample of solar-analog 
stars. This project was originally started for the purpose 
of clarifying the behavior of Li abundances (A(Li)), and 
revealed that the rotational velocity (w c sini) is the most 
influential key parameter (Takeda et al. 2007; hereinafter 
referred to as Paper I). In a successive study, we inves- 
tigated the activities of these solar analogs by using the 
residual flux at the core of the Ca II 8542 line (r (8542)), 
and confirmed a clear correlation between A(Li), w e sini, 
and ro (8542), as expected (Takeda et al. 2010; hereinafter 
referred to as Paper II). However, it turned out hard to 
discriminate the differences in ro(8542) when the activity 
is as low as that of the Sun, since it tends to get settled 
at ^0.2 and does not serve as a sensitive indicator any 
more. Actually, test calculations of non-LTE line forma- 
tion suggested (cf. Appendix B in Paper II) that the core 
flux of the Ca II 8498/8542/8662 triplet lines are rather 
inert to the chromosphcric temperature rise in the up- 
per atmosphere when the activity is low, but that for the 
Ca II 3934/3968 doublet (H+K lines) is still sensitive and 
thus more advantageous for detecting the low-level activ- 
ity. Then, we recently studied the Be abundances of these 
sample stars by using the Be II 3131 line based on the 



near-UV spectra obtained with Subaru/HDS (Takeda et 
al. 2011; hereinafter referred to as Paper III). Since these 
HDS spectra fortunately cover the Ca II H+K lines in the 
violet region, we decided to reinvestigate the activities of 
these 118 Sun-like stars by measuring the core-emission 
strength of the Ca n 3934 (K) line, in order to clarify the 
activity status of our Sun in comparison with similar so- 
lar analogs, where special attention was given to correctly 
removing the background line profile computed from the 
solar photospheric model, while taking advantage of the 
fact that atmospheric parameters of all these targets are 
well established. The purpose of this paper is to report 
the outcome of this investigation. 

2. Basic Observational Data 

2.1. Target Sample 

We use the same targets (118 solar analogs) as used in 
Papers I— III, which were selected by the criteria of hav- 
ing B — V and My values sufficiently similar to those of 
the Sun (\A(B-V)\ < 0.2-0.3 and \AM V \ < 0.3). See 
section 2 in Paper I for a detailed description about the 
sample selection. We also determined the atmospheric 
parameters (T e g, log g, vt, and [Fe/H]) from the equiv- 
alent widths of Fe lines, and the stellar parameters (M 
and age) by comparing the positions on the HR diagram 
with the theoretical evolutionary tracks (cf. section 3 
in Paper I). These parameters for most of the targets 
were actually confirmed to proximally distribute around 
the solar values (i.e., |AT cff | < 100 K, |Alogs| < 0.1 dex, 
|Aw t | <0.2 km s-\ | A [Fe/H] | < 0.2 dex, |AM| <O.1M 01 
and | A log age | < 0.5 dex; cf. figures 4 and 5 in Paper I). 

However, the following characteristics regarding the re- 
lations between these parameters are to be noted, which 
we had better bear in mind in discussing the behavior of 
stellar activities. 

— Since the effect of a decreased metallicity on B — V is 
compensated by a lowering of T g, several outlier stars 
with appreciably lower T e g as well as [Fc/H] (AT e g- < 
—200 K and [Fc/H] < —0.4 dex) are included in our sample 
(cf. figure 4c in Paper I), such as HIP 26381, HIP 39506, 
HIP 40118, and HIP 113989; and they belong to the old- 
est group (age ~ 10 10 yr). 

— These parameters are not independent from each other 
and some correlations appear to exist between specific 
combinations; such as age vs. T e g (lower T e g stars tend 
to be older), age vs. [Fe/H] (lower [Fe/H] stars tend be 
older), M vs. age (lower-mass stars tend to be older), and 
Vt vs. T c g (i> t tends to decrease with a lowered T e g), as 
recognized from figure 5 or figure 10 in Paper I, though 
the existence of outlier stars mentioned above partly plays 
a role in these tendencies. 

It should also be remarked that the data for HIP 41484 
given in Paper I were incorrect, as reported in appendix 
A of Paper II, where the correct results derived from a 
reanalysis of this star are presented. 
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2.2. Observational Material 

The spectroscopic observations of these 118 solar 
analogs and Vesta (substitute for the Sun) were carried 
out on 2009 August 6, 2009 November 27, 2010 February 
4, and 2010 May 24 (Hawaii Standard Time), with the 
High Dispersion Spectrograph (HDS; Noguchi et al. 2002) 
placed at the Nasmyth platform of the 8.2-m Subaru 
Telescope atop Mauna Kea, by which we obtained high- 
dispersion spectra covering ~ 3000-4600 A with a resolv- 
ing power of R ~ 60000. See section 2 of Paper III and 
electronic table El therein for details of the observations 
and the data reduction, as well as the basic data of the 
spectra (e.g., observing date, exposure time, S/N ratio at 
A~313lA). 

The counts of raw echelle spectra (including the effect 
of blaze function) around A ~ 3950A (the broad maximum 
between the two large depressions of Ca n 3934 (K) and 
3968 (H) lines) are typically about ~ 15 times as large as 
those at the UV region of A ~ 313lA, while the counts 
at core of the Ca n K line at 3934 A is about ~ 10% of 
those at A ~ 3950 A. Therefore, the S/N ratio at the deep 
absorption core of the K line is not much different from 
the value at A ^ 3131 A given in electronic table El of 
Paper III; that is, on the order of S/N ~ 100. 

3. Core Emission Measurement 

3.1. K line of ionized calcium 

Most activity studies of solar-type stars so far based on 
the core emission strengths of Ca n resonance lines ap- 
pear to utilize both K (3934 A) and H (3968 A) lines, pre- 
sumably due to the intention of reducing the systematic 
errors by averaging both two, since measurements tend to 
be done rather roughly by directly integrating raw spectra 
at the specified wavelength regions. 

In this investigation, however, we focus only on the for- 
mer K line at 3933.66 A, since (1) it is by two times 
stronger than the H line and thus comparatively more 
suitable as a probe of the condition at the optically- 
thin chromospheric layer, and (2) the latter Ca II H line 
at 3968.47 A is blended with the Balmer line (He at 
3970.07 A) which would make the situation more com- 
plicated in simulating the photosphcric profile to be sub- 
tracted. 

Our spectra in the selected wavelength region (3930- 
3937 A) including the relevant Ca n K line are shown 
in figures 1 (Vesta/Sun), 2 (all 118 stars including 
Vesta/Sun), and 3 (all stars overplotted) , where the con- 
tinuum levels of the observed spectra are so adjusted as 
to match the theoretical ones as explained below. We can 
see from these figures that the strengths of core emission 
considerably vary from star to star, while that for the Sun 
is apparently weak. 

3.2. Photo spheric "profile matching 

Given that the very strong Ca n K and H lines with 
considerably extended damping wings are predominant 
at ~ 3900-4000 A, it is hopeless to empirically cstab- 
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Fig. 1. Chromospheric emission feature at the core of the 
Ca II 3933.663 line for the case of the Sun (Vesta), in compar- 
ison with the theoretical photosphcric spectra computed from 
the classical model atmosphere. Three solid lines indicate the 
theoretical residual flux spectra (normalized by the contin- 
uum flux), r-^ h (= /'',x,,,t); which are based on essentially 
the same Kurucz's (1993) ATLAS9 solar atmospheric model 
but with different surface locations (i.e., the optical depth 
at the first mesh point), as depicted in the T(t5ooo) struc- 
tures shown in the inset. Red line is the LTE line profile 
derived from the adopted model atmosphere with its surface 
at logT5ooo = —5, while blue lines are the LTE and NLTE line 
profiles derived from the specially-extrapolated model atmo- 
sphere with its surface at logrsooo = —7- Open symbols ■■■ 
observed spectrum of Vesta/Sun, r^ bs (= D^ bs /ZJ°^ t ), where 
f)° bs is the actually recorded spectrum (ADU counts) while 
D™ t (regarded as an adjustable free parameter) was ad- 
equately chosen by requiring a satisfactory match between 
r^ 1 and r^ bs within < 2-3 A from the line center (except- 
ing the core-emission region) as explained in subsection 3.2. 
The specified line-core region (3932.8-3934.6 A) is indicated 
by a horizontal bar, over which the integration was made for 
evaluating the total chromospheric emission. 




theoretical photospheric 
spectra (offset by -0.05) 
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Fig. 3. In blue lines are overplotted the observed spectra 
(r° bs ) of Ca II 3933.663 in the 3930-3937 A region for all 118 
stars (the same as those shown in figure 2), along with the 
Vesta/Sun spectrum highlighted in the red thick line. The 
theoretical photosphcric spectra (r^ h ) for each of the stars 
are also overplotted in green lines (only the solar spectrum 
is in red) with a downward offset by —0.05. The wavelength 
scale of all spectra is adjusted to the laboratory frame. 
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lish the continuum position from the observed spectrum 
£)° bs (where an cchelle order covers only ^50 A). We 
thus "adjusted" the continuum position (-D°ont) °f * ne 
spectrum (judged by eye-inspection) in such a way that 
r° x hs (= D° x hs /D°^ t ) satisfactorily matches the theoreti- 
cally calculated residual flux r\ h (= F^V-^cont) 1 m tne m " 
ner wms of the line (within |AA| < 2-3 A fr om the line 
center, excepting the core-emission region). An example 
of such an accomplished match is displayed in figure 1 for 
the case of the Sun (Vesta) . 

Regarding the computation of theoretical spectra, we 
used Kurucz's (1993) WIDTH9 program, which was mod- 
ified by Y. Takeda to enable spectrum synthesis by in- 
cluding many lines. As to the atomic line data, we in- 
voked Kurucz and Bell's (1995) compilation and included 
all available lines in the relevant region. In particular, the 
data for the Ca II K line at 3933.663 A of our primary 
concern are as follows; xiow = 

0.00 cV, logs/ = +0.134, 
logTR = 8.20 [radiation damping width (s" 1 )], logr e /iV e = 
—5.52 [Stark effect damping width (s" 1 ) per electron den- 
sity (cm" 3 ) at 10 4 K], and \ogT w /N R = -7.80 [van der 
Waals damping width (s _1 ) per hydrogen density (cm -3 ) 
at 10 4 K]. Since the atmospheric parameters (T c g, logs, 
[Fe/H], and v t ) are already established (as summarized 
in table 1), the model atmosphere for each star was gen- 
erated by 3-dimensionally interpolating Kurucz's (1993) 
ATLAS9 model grids (LTE, plane-parallel model) in terms 
of T e g, logs, an d [Fe/H]. Then, the synthetic spectrum 
was computed by using the relevant atmospheric model 
along with the metallicity-scaled abundances (for all el- 
ements including Ca; i.e., [X/H] = [Fc/H] for any X) as 
well as the microturbulence (v t ), and further broadened 
according to the macrobroadening parameter determined 
in Paper III. 

3.3. Evaluation of chromo spheric emission 

Now that the theoretical photospheric background pro- 
file (r\ h ) used for subtraction has been successfully fitted 
with r^ bs by appropriately adjusting the continuum po- 
sition, we can calculate the absolute emission flux (F^ p ) 
at the K-line core originating from the chromosphere (i.e., 
after subtraction of the photospheric component) as 

Fk P = Fott[ \r? s -r i £)d\ (1) 

where Ai and A2 defining the integration range were cho- 
sen to be 3932.8 and 3934.6 A, respectively (cf. figure 
1). In order to demonstrate how this subtraction process 
works well, the photospheric profile (r\ h ) at [Ai, A2] is dis- 
played by red dashed line (along with the observed spec- 
trum shown by symbols) for each star in figure 2. Finally, 
we can obtain R' Kp (the ratio of the chromospheric emis- 
sion flux at the K line to the total bolometric flux) as 

1 We here use the astrophysical flux (F) defined by i~F\ h = 
2tt J fil^ h (0, n)dfi, with which the effective temperature T c g is 
related as 7r_Fbol = t J q F\d\ = crT^ (a: Stephan— Boltzmann 
constant). 



R' Kp ^F^/F hol = TrF^/(aT c 4 s ). (2) 

The resulting log-R^p values for each of the 118 solar 
analogs (+Sun) are presented in table 1, where other 
activity-related quantities (ro(8542) and w sini) deter- 
mined in Paper II are also given, along with the Li/Bc 
abundances and stellar parameters established in Papers 
I and III. 

3.4- Zero-Point Uncertainties in R' Kp 

We would like to remark here that the "absolute" values 
of R' Kp are not very meaningful in the low-activity regime 
because of the uncertainties in its zero-point. That is, the 
theoretical profile we have computed for subtraction of 
the background photospheric component is by no means 
uniquely defined. Actually, apparently different results 
may be obtained depending on how it is calculated. 

This situation is demonstrated in figure 1 for the solar 
case. The solar atmospheric model (which was obtained 
by interpolating the grids of Kurucz's ATLAS9 model at- 
mospheres) we adopted for calculating the LTE photo- 
spheric profile (red line) has its surface at log T5000 = ~5 
(T 8urf ~ 4000 K). However, we note that the same but 
simply extrapolated model up to logrfQQQ = — 7 (T surf ~ 
3600 K) yields an appreciably deeper core (upper blue 
line), reflecting the fact that the residual flux at the cen- 
ter is determined by - B x (T surt )/B x (T ph ) (T ph : photo- 
spheric temperature) and very sensitive to T s in this 
violet region where Wicn's approximation nearly holds. 
Moreover, the core of the non-LTE profile (simulated by 
using the departure coefficients computed in Paper II for 
the case of Model E; cf. Appendix B therein) gets even 
more deeper approaching a completely dark core (lower 
blue line). Besides, the core shape strongly depends on 
the turbulent velocity field in the upper atmosphere; e.g., 
compare the non-LTE profile in figure 1 (computed with 
a depth-independent microturbulence of 1 km s -1 ) with 
that for Model E depicted in figure B.l(b) of Paper II 
(where a variable microturbulent velocity field increasing 
with height was adopted). 

Accordingly, R' Kp values [equation (2)] may be un- 
certain by an arbitrary constant, because r\ h d\ in 
equation (1) can have different values depending on how 
r\ h is computed, though "relative" differences between 
^kp values of each star are surely meaningful as long as 
they are evaluated in the same system. Therefore, we 
should keep in mind that comparison of absolute values 
of our R' Kp with those of similar R' parameters derived 
by other groups, which we will try in subsection 4.1, is 
not much meaningful when the low-activity region (e.g., 
logR^p < —5) is concerned. (On the other hand, such 
a zero-point problem should not be so serious for higher- 
activity cases, where logi?^ p tends to be dominated by the 
emission component and the role of r x h subtraction is less 
significant.) In any event, we consider that our choice of 
shallower r\ h is adequate, since it leads to smaller values 
of i?^p (as a result of larger subtraction), which even- 
tually realizes a larger contrast in the logi?^ values of 
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low-activity stars. 

4. Discussion 

Comparison of R' Kp with previous studies 

The logR^p values determined in subsection 3.3 are 
compared with the equivalent activity indices 2 derived by 
Strassmeier et al. (2000) [logi?y, Wright et al. (2004) 
[logi?jj K ], and Isaacson and Fischer (2010) [logi?y K ] in 
figures 4a, 4b, and 4c, respectively. These figures re- 
veal almost the same tendency of correlation between 
our logi?' and those of three studies: The agreement 
at logi?' > —5 is mostly good, though our logi?' tends 
to be slightly larger by ~ 0.1 dex for high-activity stars 
of logi?' > —4.5. Meanwhile, a distinct difference is ob- 
served at the low-activity region where their logi?' values 
tend to settle down at ~ —5, while ours are dispersed 
over —5.5 < logi?' < —5. Admittedly, there is not much 
meaning in comparing the absolute i?' values of different 
systems with each other, as remarked in subsection 3.4. 
The important point is, however, that we could detect the 
subtle difference in the low-level activity by measuring the 
weak core-emission at Ca II K with a careful subtraction 
of the background photospheric profile, while such a pre- 
cision for distinguishing the delicate difference in the weak 
emission strength could not be accomplished by compar- 
atively rough measurements in those previous studies. 

4-2. Connection with stellar parameters 

Figures 5a, 5b, and 5c display how the three activity- 
related parameters (ro(8542), u c sini, and log age; cf. 
Paper II) are correlated with logi?^ p . 

We can observe in figure 5a a marked sensitivity- 
difference between log i?^ p and ro(8542) in the low- 
activity region; i.e., the latter is inert to a variation of 
low-level activity and stabilizes at ~ 0.2, despite that the 
former still shows an appreciable variability over — 5.5 ^ 



logi?^ p < — 5. This is just what we have expected (cf. 
Appendix B in Paper II), and demonstrates the superi- 
ority of the Ca II HK core emission (as long as correctly 
measured) to the line-center residual flux of Ca n 8542 
when it comes to investigating the activities of solar-type 
stars as low as the Sun. 

Figure 5b shows a positive correlation between R' Kp and 
w sini, suggesting that stellar activity depends on the ro- 
tation rate, as we already confirmed in Paper II (cf. figure 



Since Strassmeier et al. (2000) treated H and K lines sepa- 
rately and presented each data of and F^r , we could convert 
their i?Jj K [= (FL + FL)/Fbol] ( as clearly defined by them) into 
R' K (= F^/Fbal) which is directly comparable with our R' Kp - 
Meanwhile, the R'n K data derived by Wright et al.'s (2004) as 
well as Isaacson and Fischer (2010) appear to be essentially the 
average (not the sum) of R' n and R' K as judged by their extents. 
Therefore, we should be cautious about different definitions in 
the meaning of R'^j^- This situation is manifestly displayed in 
figure 7 of Paper II, where we can see that R'n K (Strassmeier) is 
systematically larger than R'n K (Wright) by 0.3 dex. Therefore, 
we compared our log-R^ p with /?^ K (Wright, Isaacson) without 
applying any correction, since the latter is practically equivalent 
to R' K in any case. 
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Fig. 4. Correlation of the logiJ^, p indices determined in this 
study with the literature values taken from three represen- 
tative papers: (a) Strassmeier et al.'s (2000) logR' K values, 
where 16 stars are in common, (b) Wright et al.'s (2004) 
logiJJj K values, where 50 stars are in common, (c) Isaacson 
and Fischer's (2010) logi?Jj K values, where 66 stars are in 
common. (Since logiiJj K of Wright et al. (2004) as well as 
Isaacson and Fischer (2010) should be equivalent to the mean 
of log R' u and log R' K , it may be directly compared with our 
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Fig. 5. Diagrams showing how the activity-related quanti- 
ties derived in Paper I and Paper II (ro(8542), v c sini, and 
log age; also given in table 1) are correlated with the activ- 
ity index (logi?^) derived in this study, (a) ro(8542) vs. 
logJ?^ p , (b) t> c sini vs. logi?^ p , and (c) log age vs. log_R^ p . 
The log R' Kp q value which we derived from the spectrum 
of Vesta/Sun (corresponding to the near-minimum phase of 
activity) is indicated by the bigger (red) circle, while the 
expected minimum-maximum range of log R' Kp q (between 
—5.35 and —5.15; cf. subsection 4.3) is shown by a horizontal 
bar. 



5a therein). However, since v e sini values cluster around 
~ 2 km s _1 at — 5.5 ^ logi?^ p < — 5, we can not state much 
about whether this activity-rotation connection persists 
down to such a low-activity region (also, uncertainties in 
the projection factor prevent from a meaningful discus- 
sion) . 

When we compare the age vs. i?^ p relation depicted 
in figure 5c with the similar ro(8542) vs. age plot (cf. 
figure 5c of Paper II) , the anti-correlation is more clearly 
(or less unambiguously, to say the least) recognized in the 
present case, thanks to the extended dynamic range of 
the activity indicator for low-activity stars, though the 
dispersion is still considerably large. 

How the abundances of Li and Be depend on R' Kp de- 
termined in this study is illustrated in figure 6, where 
their dependences upon r (8542) and v e suii already dis- 
cussed in Papers II and III are also shown for comparison. 
We can see from figures 6a that the near-linear relation 
between A(Li) and logi^p (^( Li ) - 7 + l°g- R K P ) holds 
widely from highly active (logi?^ p ~ —4) to less active 
(log-Rj^p ~ —5.5) stars, in contrast to the case of ^4(Li) vs. 
ro(8542) where ^4(Li) shows a considerable dispersion at 
r (8542) ~ 0.2 (as if compressed) because of the less sen- 
sitivity of ro (8542). This substantiates the observational 
conclusion in Paper II (or corroborates its validity even 
for less-active cases) that ^4(Li) closely depends upon the 
stellar activity, which further lends support for our previ- 
ous argument that the key parameter for controlling the 
surface lithium in solar-analog stars is the stellar rotation. 

Regarding A(Be), we suspected in Paper III that the pe- 
culiar 4 stars showing drastically depleted Be (by > 2 dex) 
in comparison with the solar abundance (while others have 
more or less near-solar Be) might be very slowly-rotating 
and and thus less active star than the Sun. However, since 
these 4 stars have appreciably different logi?^ p from each 
other (—5.4 < logi?^ p < —4.9; cf. figure 6a') and are not 
necessarily less active than the Sun, this speculation docs 
not seem likely. Some other explanation would have to be 
sought for. 

4-3. Activities of solar analogs and the Sun 

We now discuss the activity trends of 118 solar analogs, 
with special attention paid to the status of our Sun among 
its close associates, as the main subject of this study. 
While our discussion is based on the log-R^p indices de- 
termined by ourselves, we should keep in mind that only 
one snapshot data is available for each star. This means 
that uncertainties due to possible time- variations of stel- 
lar activities (whichever cyclic or irregular) are inevitably 
involved. For example, the Mt. Wilson S value for the 
Sun varies over the range of 0.16 < Sq < 0.20 (Baliunas 
et al. 1995), which may be translated into a variability in 
logi?^ by - 0.2 dex (-5.0 < logi?^ < -4.8) according to 
the transformation formula given by Noyes et al. (1984). 
Since the observation time (2010 February 5) of our solar 
spectrum (Vesta) had better be regarded as correspond- 
ing to the near-minimum phase because of the apprecia- 
bly retarded beginning of cycle 24 after the minimum in 
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Fig. 6. Abundances of Li (left panels) and Be (right panels) 
plotted against logi?^ (top), ro(8542) (middle), and v c s\ni 
(bottom). The solar values are indicated by the bigger (red) 
circle in each panel. The expected minimum-maximum range 
of logfj^p g is shown by a horizontal bar. The results for the 
determinable cases are shown by filled circles, while the open 
inverse triangles denote the upper limits for the unmcasurablc 
cases. 

2008, our \ogR' Kp Q (= —5.33) may as well be raised by 
^$0.2 dex at the solar maximum phase. It is thus rea- 
sonable to assume that our log-R^p Q ranges from —5.35 
(solar minimum) to —5.15 (solar maximum), as indicated 
by a short horizontal bar in figures 5,6, and 7. 

The distribution histogram of logR^ for all the 118 
stars (and the Sun) and a similar histogram of ro(8542) 
(for comparison; cf. Paper II) arc shown in figures 7a 
and 7b, respectively We immediately notice in figure 7a 
the bimodal distribution of log R' Kp having two peaks at 
~ —5.3 and ~ —4.3, constituting a well-known Vaughan- 
Preston gap (Vaughan & Preston 1980), while this bi- 
modal trend is not clear in the distribution of ro(8542) 
(figure 7b) because of the densely peaked population 
around ro(8542) ~ 0.2, reflecting its insensitivity when 
the activity is low (cf. figure 5a). 

As we can see from figure 7a, the Sun with logi?J {p Q of 
—5.33 manifestly belongs to the low-activity group (rang- 
ing from ~ —5.4 to ~ —5.0). As a matter of fact, only 
11 (~ 10%) 3 out of 118 solar analogs have logi?j <p values 



This ratio naturally depends on the reference solar activity, 
which may be subject to uncertainties due to cyclic variations 
as mentioned at the beginning of subsection 4.3. For example, if 
we assume a somewhat higher value of —5.2 for logi?^ p g cor- 
responding to the active phase of the Sun, this fraction becomes 
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Fig. 7. Histograms showing the distributions of (a) logi?^, p 
and (b) ro(8542) for our sample of 118 solar analogs. The 
position for the solar value is indicated by an downward arrow 
at each panel. The expected minimum-maximum range of 
logi?^, p q is shown by a horizontal bar. 

smaller than —5.33. This distinctly low-activity nature of 
the Sun is also recognized by eye-inspection of figure 2, 
revealing that the Ca n K line emission strength in our 
solar spectrum is near to the minimum level among other 
stars. 

Thus, we can rule out the possibility for the existence 
of a significant fraction of Maunder- minimum stars (i.e., 
solar-type stars with appreciably lower activity than the 
Sun, even showing an another peak well below the cur- 
rent solar- minimum level), such as those once suggested 
by Baliunas and Jastrow (1990). This result corrobo- 
rates the arguments raised by recent studies (e.g., Hall & 
Lockwood 2004; Wright 2004), which cast doubts about 
the reality of such a high frequency of Maunder-minimum 
stars. Thus, our Sun belongs to the group of manifestly 
low activity level among solar analogs, the fraction of stars 
below which is essentially insignificant. 

4-4- Stars of Subsolar Activity 

Although we have concluded that the Sun belongs to 
nearly the lowest activity group, some stars do exist show- 
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ing activities still lower than that of the minimum-Sun, 
which are worth being examined more in detail. Since 
we defined that solar logi?^ p Q varies from —5.35 (mini- 
mum) to —5.15 (maximum) (cf. section 4.1), such stars 
may be sorted out by the criterion of log-R^p < —5.35, 
which resulted in the following 8 objects (logi?^ p , A(Li), 
log age, [Fe/H], remark): HIP 7918 (-5.42, 1.89, 9.56, 
+0.01), HIP 31965 (-5.39, 1.00, 9.87, +0.05), HIP 39506 
(-5.44, < 0.9, 10.36, -0.62) HIP 53721 (-5.43, 1.75, 
9.89, -0.02, PHS), HIP 59610 (-5.40, 1.62, 9.63, -0.06, 
PHS), HIP 64150 (-5.38, < 1.0, 9.63, +0.05, Be depleted), 
HIP 64747 (-5.40, < 1.1, 9.80, -0.18), and HIP 96901 
(-5.38, < 1.1, 9.77, +0.08, PHS). 

While only one (HIP 39506) of these is an outlier of 
lower [Fe/H] as well as lower T e s (belonging to rather old 
population) as mentioned in subsection 2.1, the remaining 
7 stars arc Sun-like stars with sufficiently similar param- 
eters, which excludes the possibility of such a low-level 
activity being due to stellar-evolution (i.e., evolved sub- 
giants; cf. section 1). 

We note here that (1) three (out of only five in our 
sample of 118 stars) planet-host stars are included, (2) all 
these stars have low-scale Li abundances (A(Li) < 2), and 
(3) their ages are similar to or older than that of the Sun. 
Combining these facts with the consequences in Papers I 
and II, we consider that these stars have actually low ac- 
tivities even compared with the solar-minimum level, and 
this is presumably attributed to their intrinsically slow ro- 
tation (which is closely related with the Li abundance as 
well as with the existence of giant planets). 

It is, therefore, interesting to investigate the variabil- 
ities of these low-activity stars by long-term monitoring 
observations, in order to see how their activities behave 
with time (cyclic? flat? irregular?). Admittedly, activity 
observations for these stars have been reported in sev- 
eral published studies so far, for example: HIP 7918, 
53721, 64150, and 96901 by Duncan et al. (1991); 
HIP 53721, 59610, 64150, and 96901 by Wright et al. 
(2004); HIP 7918, 53721, and 96901 by Hall et al. (2007); 
HIP 31965, 59610, 64150, and 96901 by Isaacson and 
Fischer (2010). However, as they are still quantitatively 
insufficient for establishing the long-term behavior of their 
activities , much more observations are evidently needed. 

5. Conclusion 

There have been several arguments regarding the status 
of solar activity among similar Sun-like stars, which be- 
gan with the implication of Baliunas and Jastrow (1990) 
based on their Mt. Wilson HK survey project that a con- 
siderable portion (~ 1 /3) of solar-type stars have activities 
significantly lower than the present-day Sun, which they 
called "Maunder-minimum stars." However, their con- 
clusion could not be confirmed by Hall and Lockwood's 

4 As a comparatively well studied case where a sufficient amount 
of data are available, we may presumably state that HIP 96901 
= HD 186427 showed a Sun-like cyclic variation in the 1995- 
2007 period; cf. figure 7 of Hall et al. (2007). 



(2004) follow-up study, and Wright (2004) criticized the 
reality of such considerably low-active solar-type stars by 
pointing out that most of them are not so much dwarfs as 
evolved subgiants. 

Given this controversial situation, we decided to con- 
tend with this problem by ourselves based on carefully 
selected sample of 118 solar-analogs sufficiently similar to 
each other (which we already investigated their stellar pa- 
rameters as well as Li/Be abundances in a series of our 
previous papers), with a special attention being paid to 
reliably evaluating their activities down to a considerably 
low level. 

Practically, we measured the emission strength at the 
core of Ca II 3933.663 line (K line) on the high-dispersion 
spectrogram obtained by Subaru/HDS, where we gave ef- 
fort to correctly evaluating the pure emission component 
by removing the wing-fitted photospheric profile calcu- 
lated from the classical solar model atmosphere, which 
enabled us to detect low-level activities down to logi?^ p ~ 
-5.5. 

A comparison of our log R' Kp results with the corre- 
sponding logi?' values of Strassmeier et al. (2000), Wright 
et al. (2004), and Isaacson and Fischer (2010) revealed 
that low-active stars (for which they derived logJ?' ~ —5.1 
at the minimum limit) actually have a dispersion of ~ 
0.4 dex (—5.5 < logi?^ p < —5.0) in our measurement, 
suggesting that our log-R^p has a higher sensitivity and 
thus advantageous. A similar situation holds regarding 
the comparison with ro(8542) we used in Paper II; i,e., 
this index stabilizes at ~ 0.2 and becomes insensitive for 
low- active stars in contrast to logi?^ p . 

As another merit of using logi?^ p , we can state that the 
visibility of the A(Li)-activity relation as well as the age- 
activity relation becomes comparatively clearer, because 
this activity index turns out to have well diversified values 
for low-activity stars thanks to its high sensitivity, which 
can not be accomplished by using, e.g., ro(8542). 

From the distribution histogram of logi?^ p , we could 
recognize a clear Vaughan-Preston gap between two peaks 
at ~ —5.3 and ~ —4.3. Our result of logi?^ Q = —5.33 
manifestly suggests that the Sun belongs to the group of 
the former peak and has a distinctly low-active nature 
among solar analogs. Actually, a fraction of stars with 

log-R^p — l°g^kp.o 18 onrv ~ 10%. This consequence ex- 
clude the possibility for the existence of a considerable 
fraction (e.g., ~ 1/3) of "Maunder-minimum stars" such 
that having activities significantly lower than the current 
solar-minimum level as once suggested by Baliunas and 
Jastrow (1990). 

Yet, some stars (only a minor fraction of the sample) do 
exist showing activities still lower than that of the solar- 
minimum level. Having examined such 8 low-activity 
stars, wc found that they tend to include planet-host stars 
and have low Li abundances, from which we suspect that 
their activities are actually low as a result of intrinsically 
slow rotation. It would be an important task to clarify the 
behavior of their activity variations by long-term monitor- 
ing observations. 
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Fig. 2. Display of the 3930-3937 A region spectra (r^ bs ) of the Ca II K line at 3933.663 A for all the 118 program stars along 
with the Vesta/Sun, where the theoretical line profile calculated from the solar photospheric model is also overplotted by (red) 

dashed line in the core region (3932.8-3934.6 A) around the line center. Each spectrum is vertically shifted by 0.1 (in continuum 
unit) relative to the adjacent one. The wavelength scale of all stellar spectra is adjusted to the laboratory frame by correcting the 
radial velocity shifts. The HIP numbers are indicated in the figure. 
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Table 1. Activity index, rotation, Li abundance, age, and the atmospheric parameters. 
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Table 1. (Continued) 
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Note. Following the HIP number in column 1, logR' Kp (column 2) is the emission-strength index of the Ca II K line determined in this 
study, ro(8542) (column 3) is the residual flux at the line center of the Ca II 8542 line (cf. Paper II), i; e sini (column 4) is the projected 
rotational velocity in km s _1 calibrated in Paper II, A(Li) (column 5) as well as A(Be) (column 6) are the lithium and beryllium abundances 
derived in Papers I and III (in the usual normalization of A(H) = 12.00), and log age (column 7) is the logarithm of the stellar age (in yr) 
derived from theoretical evolutionary tracks in Paper I. In columns 8-11 are presented the atmospheric parameters for each star, which arc 
the "standard parameters" derived from Fe I and Fe II lines (cf. section 3.1.1 in Paper I): T e ff (effective temperature: in K) logcj (surface 



gravity; in cm s ), Vt (microturbulence; in km s 
abundance), respectively. In column 12 (remark), 



1 ), and [Fe/H] (metallicity defined by A(Fe)- 
'PHS" denotes planet-host stars. 



7.50, where A(Fe) is the logarithmic Fe 



